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Introduction 


Finite-element models are used extensively by the rotorcraft industry for determining the vibratory 
characteristics of helicopter airframes. These airframe models are typically large, usually containing several 
thousand finite elements, and their definition involves considerable effort. Results obtained by using these 
models usually do not correlate well with experimental data, particularly at the higher frequencies that are 
near the predominant rotor excitation frequencies (refs. 1 through 4). Although there are several possible 
reasons for the lack of correlation, one that has recently attracted attention is the accuracy of the elements 
themselves. 

Most finite-element codes in use today employ what are commonly referred to as h-version elements. An 
h-version element uses fixed-order shape functions to relate the discrete nodal displacements of the element 
to the continuous displacements within the element. Mesh size normally controls the accuracy of the results 
obtained using h-version elements. As the number of elements is increased, the mesh size decreases and the 
results tend to converge. Tests for convergence usually involve solving a series of problems with successively 
refined meshes until convergence is indicated. Such convergence checks are often used in studies which employ 
small models. To check the finite-element model of a large structure (such as an airframe) for convergence by 
mesh refinement would be a formidable task not likely to be done in practice. So-called p- version elements that 
have recently gained attention have the potential for both improving accuracy and simplifying convergence 
checks. 

The p-version element is different from the h-version element in that the user may select the order of the 
shape function defining the displacement behavior of the element. In general, the higher the order of the shape 
function, the more accurate the results. Increasing the order of the shape functions is analogous to increasing 
the number of elements for a comparable h-version model. This feature becomes a distinct advantage in 
convergence checks of large finite-element models because raising the order of the shape functions would be 
much easier than increasing the number of elements in the data input file. 

The advantage of using p-version elements is lost if the element does not closely approximate the geometry 
of the structure being modeled. For example, a uniform p-version beam element would pose no advantage 
over a uniform h-version element in either accuracy or ease of use in modeling a tapered beam. However, 
a p-version beam element whose cross-sectional dimensions vary linearly with length would allow for a more 
accurate representation of a tapered beam than could be achieved with a uniform beam element. Because 
many substructures in a typical airframe consist of beams whose cross sections vary along their length, the use 
of p-version elements offers an advantage over h-version elements only if the p-version elements can accurately 
model the physical structure without resorting to excessive mesh refinement. 

The Langley Research Center has underway a research structural dynamics program which is investigating 
finite-element modeling techniques for helicopter vibrations analyses. This activity has identified the need 
for a small research code containing a limited library of p-version beam and shell elements. To meet the 
requirements of a research tool, these elements must be well documented, easily modified, and tailored to 
the needs of airframe vibrations analyses. Work has been initiated at Langley to develop such a code. Initial 
attention is being given to the beam element. A survey of the finite-element literature revealed no commercially 
available codes with p-version beam elements and only a few research codes having such elements. None of the 
latter codes appear to be suitable for airframe vibrations analyses work. However, there has been significant 
work done in the area of tapered h-version beam elements which is relevant to the beam element that is the 
subject of this paper. Pertinent work in this area is summarized. 

Tapered beam finite elements based on cubic polynomial shape functions have been under investigation 
for many years (refs. 5 and 6). Higher order, tapered beam elements have also been developed. Thomas and 
Dokumaci (ref. 7) presented two tapered beam elements using quintic polynomials. One beam element had 
two internal nodes, whereas the other beam element defined the second derivative at each end of the beam as 
a degree of freedom. Interelement continuity was then enforced for displacement and the first two derivatives, 
which formed a C 2 -type element. A tapered beam element using seventh-degree polynomial shape functions 
was derived by To in reference 8. Each node had four degrees of freedom (displacement and the first three 
derivatives) on which interelement continuity was imposed. It should be pointed out that these higher order 
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elements are not p~ version elements because the order of the shape functions was fixed when the element was 
derived. In another paper, To (ref. 9) developed a tapered beam element incorporating shear deformation and 
rotary inertia but based the element formulation on cubic shape functions. More recently, a tapered p-version 
beam element was developed by Hodges, Hopkins, Kunz, and Hinnant (ref. 10) specifically for modeling rotor 
blades. To correctly model a spinning rotor blade, the element includes the nonlinear effects of large nodal 
displacements and rotations, aerodynamics, and inertial rotation. Such a complex element is not needed for 
modeling airframe structures. As a practical matter, elements used for the analysis of airframes should be 
relatively simple because of the large number of elements normally required to model the structure adequately. 

The objective of this paper is to present the derivation of a tapered, p-version beam element for use in 
dynamic analyses of general structural systems. This element features hierarchical shape functions which allow 
higher order analyses to use element matrices established for lower order analyses. Appropriate orthogonal 
relations for the shape functions are employed to avoid ill-conditioned matrices and reduce the number of 
nonzero terms. Element matrices are explicitly formed, thus eliminating the need for numerical quadrature 
resulting in a simpler implementation and a reduction in roundoff error. 


Symbols 

A 

o>i 

E 

fuiifviifoi 

G 

Idof 

IxxJyyJzz 

J 

K 

Kee 

x ru tsv \sw ttO 
^EI^EI^EI^EI 

\SU TSV T {W T s0 

l 

M 

M ee 

m y,, m;,, Mf,, Mf, 

Ndof 

Nelem 

iVf, NV, N?, N f 


cross-sectional area 

coefficients of polynomial representing cross-sectional property 
modulus of elasticity 

ith axial, lateral bending, and torsional frequency, 
respectively, Hz 

shear modulus 

number of internal degrees of freedom 

area moment of inertia about X-, T-, and Z-axis, respectively 

torsional stiffness constant 

stiffness matrix 

external stiffness submatrix 

external-internal stiffness coupling submatrice associated with u, 
u, tu, and 0, respectively 

internal stiffness submatrice associated with u, t?, w , and 0, 
respectively 

length of beam 

mass matrix 

external mass submatrix 

external-internal mass coupling submatrice associated with ti, v, 
w, and 0, respectively 

internal mass submatrice associated with u, v, w, and 0, 
respectively 

total number of degrees of freedom 
number of elements 

ith shape function associated with u, v, w, and 0, respectively 
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7VP 

PuXvXltlXo 


P2(x) 

P*(x) 


sfS“ 

s « 

SJT 


t 

V 


u,v,w 

XXr Z 


x,y,z 

a t A 9 ’ 7? 

«i > #. l\ 
«() 


A(i) 

9 

P 


zth shape function of C° type 
ith shape function of C l type 

number of shape function associated with u, v, w , and 0, 
respectively 

quadratic polynomial representing generic cross-sectional 
property 

quartic polynomial representing generic cross-sectional property 

ith discrete degree of freedom associated with u, v, w, and 0, 
respectively 

matrix function defined by Jjj P4(x)Nf N® dx 
matrix function defined by /q P4(x)(N^Y(NjY dx 
matrix function defined by jjj P2(x)N^Nj dx 
matrix function defined by Jq P4{x)(N^) f (Njy dx 
matrix function defined by /q P4 ( x ) ( ) ff ( Nj ) 1 f dx 
kinetic energy 
time, sec 
strain energy 

displacement along X-, F-, and Z-axis, respectively 

rectangular axis system with origin at one end of beam; 

X-axis along neutral axis of beam, and F- and Z-axes oriented 
parallel to principal axes of beam 

independent variable along X-, Y-, and Z-axis, respectively 
coefficients in generating equation for Nf 
coefficients in generating equation for iV- 1 
first variation of ( ) 

fO {i^j) 

[l (i = i) 

J -1 (i = 2 or 4) 

\ 1 (Otherwise) 

rotation about X-axis (i.e., torsion) 

mass density 

ith axial, lateral bending, and torsional frequency, respectively, 
rad/sec 


Primes to a symbol denote a derivative with respect to x. A dot over a symbol denotes a derivative with 
respect to time. 
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Basic Mathematical Formulation 

The beam finite element is derived assuming that the beam behaves kinematically like a Bernoulli- Euler 
beam. Rotary inertia effects are included; however, shear flexibility is not included in the present formulation. 
The cross-sectional area can vary as a quadratic polynomial along the length of the beam, and the area moment 
of inertia can vary as a quartic polynomial. The element axis system X , Y, Z is oriented so that the X-axis is 
along the neutral axis of the undeformed beam, and the Y- and Z-axes are parallel to the principal axes of the 
cross section. (See fig. 1.) This orientation results in a zero cross product of inertia (Iyz = 0)* The continuous 
displacements u, v, u;, and 6 are assumed to be only a function of x and time. 

Hamilton’s principle is used in the derivation of this beam element. Assuming only conservative forces, 
Hamilton’s principle is stated as (ref. 11) 


[ 1 2 6{T -V)dt = 0 

Ju 


(1) 


where t\ and 1 2 represent arbitrary times at which the state of the system is known, and T and V are the 
kinetic and strain energies, respectively. The kinetic energy of the beam is found in reference 12 as 


1 = ^ j p [Ait 2 4- Av 2 + Aw 2 + Ixx ® 2 + Izz(v 9 ) 2 + Jyy(w ') 2 


dx 


and the strain energy follows similarly from reference 13 as 

1 f l 


v = \j [eA(u') 2 + EI zz (v") 2 + EI YY (w") 2 + GJ(0') 2 ] dx 


( 2 ) 


(3) 


The continuous problem associated with the continuous displacements u, v, w, and 0 is discretized by 
introducing discrete degrees of freedom q % which are related to the continuous displacements according to 


Pu 



(4) 

i=l 


Pv 


v-zm 

(5) 

i= 1 


Pu, 


s 

II 

M 

e 

-•6 

(6) 

i= 1 


Pe 


» = I 

(7) 

2 = 1 



Substituting equations (4) through (7) into equations (2) and (3) expresses the kinetic and strain energies 
in terms of discrete degrees of freedom. The application of Hamilton’s principle (eq. (1)) then leads to the 
identification of the following nonzero terms which appear in the element mass and stiffness matrices: 


Kj = f pAN?Nf dx 

M?j = P [AN? Nj + I zz (N?)'(Njy 

J = f P [aN?’N? + lYY(Nn'(N?y 
J 0 


Mr 


dx 

dx 


( 8 ) 

(9) 

( 10 ) 
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pi XX N?Nj dx 


( 11 ) 


\ 


TSU _ 

K i,j 

£ EA{N?)'{N?y dx 

(12) 

TfV _ 

[ EI ZZ {NV)"{N])" dx 

(13) 

K K = ) 

f l Elyyi Nf ) " ( Nf ) " dx 
0 3 

(14) 

Kj = 

Jl GJ(N?)'(N<>)' dx 

(15) 


Shape Functions 

The discrete degrees of freedom are divided into two sets, external and internal. The 12 external degrees 
of freedom, which are depicted in figure 2, correspond to the usual definition of the physical nodal degrees 
of freedom for a beam finite element (ref. 14). The internal degrees of freedom have no physical significance 
but are simply the coefficients of the higher order shape functions. These internal degrees of freedom serve 
to increase the accuracy of the transformation from the discrete problem having a finite number of degrees of 
freedom to the continuous problem having an infinite number of degrees of freedom. The number of internal 
degrees of freedom used in the beam element is specified by the data input file and may vary from zero to 
theoretically infinite. Setting P u = Pq = 2 and P v = P w = 4 in equations (4) through (7) will lead to the 
classical 12-degree-of- freedom beam element. Increasing any of these P’s adds internal degrees of freedom to 
the element. The shape functions N t are usually taken to be polynomials although in theory they can be any 
set of functions. 

Shape Functions for u and 9 

The shape functions N “ and Nf are identical and have C°-type continuity. That is, continuity is enforced 
across element boundaries, but continuity of the derivatives across element boundaries is not enforced. Shape 
functions satisfying <7° continuity will be denoted by Nf (i.e., Nf = Nf = Nf). The first two shape functions 
in this set are the well-known linear polynomials (ref. 15) 

*1=— 7 + 1 (16) 

^2 = 7 (17) 

The higher order <7°-type shape functions used herein were derived subject to two requirements: First, the 
<7° continuity is enforced by restricting the higher order shape functions to be zero at the element boundaries, 
and second, the set of higher order polynomials must be orthogonal with respect to their first derivative. 
Orthogonality of the first derivative was chosen over polynomial orthogonality because the element mass and 
stiffness matrices obtained by requiring first-derivative orthogonality contain fewer nonzero terms; thus explicit 
integration is facilitated. This also results in matrices which are better conditioned than those obtained from 
orthogonal shape functions. These requirements are expressed mathematically by the following equations: 


7VP( 0)=0 (i> 3) (18) 

N?(l)= 0 (*>3) (19) 

J\n?)'{NP' dx = j6ij (*> 3J>3) (20) 
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The first three higher order C°-type shape functions in the set defined by equations (18) through (20) are 


7V 3 ° = VS 
n% = VS 




n /— / X ® xV xV 

= ^ 5 F - 10 e + 6 e -7 


In general, once N 3 is known, the i + 1 shape function can be found from the recursive formula 


A? +1 = (a? + */3?)N? - 7 P/VP_ i 


u ^ ^ 


where 


0 (2i-3)v / 2T r I 

” iV 21^3 

0_ (4f-6)y / 2r^I 

Pl ilV2i^S 

o _ (i- 3)^27^! 

7i *>^*^5 


(*>3) 

(*>3) 

(i>3) 


( 21 ) 

(22) 

(23) 

/ r\ 4 \ 

(25) 

(26) 
(27) 


Shape Functions for v and w 

The shape functions N? and N™ require C 1 -type continuity, meaning that both the functions and their 
first derivatives must be continuous across element boundaries. The shape functions Nf and N™ are identical 
except that A^" = — N$ and N™ = -7V| to ensure that the discrete rotational degrees of freedom <75 and qn 
have the sense indicated in figure 2. Shape functions satisfying C 1 continuity will be denoted by AT* (i.e., 
AT = \(i)N™ = N--). The first four C^-type functions for the beam element are (ref. 15) 


Nl 

*2 

Nl 




x ^ x ^ 

J2~ 2 T + X 


X 3 x 2 

_2 ]3 +3 ]2 



(28) 

(29) 

(30) 

(31) 


The derivation of the higher order C^-type shape functions employed in this paper is based on the same 
philosophy as the derivation of the C 0 -type shape functions. To ensure C l continuity the higher order shape 
functions must have zero slope and displacement at the element boundaries. Consistent with the previous 
discussion on the orthogonality properties of the C°-type shape functions, the C^-type shape functions are 
required to be orthogonal in their second derivative. The specific requirements are 


Nh 0) = 0 


(i>5) 


(32) 
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(33) 

(34) 

(35) 


(JV/)'(0) = 0 

Nl(l) = 0 

(Nl)'(l) = 0 


{i > 5) 
{i > 5) 
{i > 5) 


J\nI)"(N})" dx = ±614 ( * > 5, j > 5) 


(36) 


The first three higher order C 1 -type shape functions in the set defined by equations (32) through (36) are 

i r- ( x 4 x 3 x 2 \ 

Nb = ^\2F~¥ + w) 


i /= / X 5 5x 4 X 3 X 2 

= ^ ( "T 3 " + W ~ 2 1 3 + 2p 


7x° 


N l = ^\w- 7 « + 


r 5 15x 4 10x 3 x 2 


+ 


/5 2/ 4 3/ 3 2/ 2 


(37) 

(38) 

(39) 


If AT 4 is given, can be found from the recursive relation 

Ni + i = (°i+*0i)*fi-1iNU (*> 5) (40) 

where 


! (2* - 5) n /2T=73 

a * iy/2i — 5 

(*>5) 

(41) 

l (4*-10)>/2t-3 
iW2i - 5 

(*>5) 

(42) 

, (i-5)v/2*-3 

ll W2i - 7 

(i> 5) 

(43) 


Matrix Functions 

Inspection of equations (2) and (3) indicates that there are eight cross-sectional properties associated with 
a beam element. Two of them, pA and EA , are represented by the quadratic polynomial 

P 2 (x) = a 2 x 2 + aix + ao ( 44 ) 

The remaining six properties, plxx * PWYi P^ZZi GJ, Elyy, and EIzz-> are represented by the quartic 
polynomial 

Pj(x) = a±x 4 + agx 3 + a 2 x 2 + aix + ao (45) 

Because these polynomials can represent any of the aforementioned cross-sectional properties by substituting 
in the appropriate values for the coefficients a^, equations (44) and (45) are referred to herein as generic 
cross-sectional property polynomials.” Note that P 2 (x) is a subset of Pi(x), and, therefore, the coefficients ao 
through a 2 refer to both polynomials. The context in which the coefficients are used will dearly specify which 
generic cross-sectional property polynomial is being referenced. 

Substituting the generic cross-sectional property polynomials for the actual cross-sectional properties in 
equations (8) through (15) and substituting generic N° or N l shape functions for the displacement-specific 
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shape functions, the 10 terms in equations (8) through (15) are reduced to 5 unique terms. These five unique 
terms are referred to herein as “matrix functions,” and have the following definitions: 


P 4 (x)N?Nf dx 

(46) 

P 4 (x)(Nf)'(N^y dx 

(47) 

P 2 (x)nIn} dx 

(48) 

P 4 (x)(NlY(N}y dx 

(49) 

P 4 (x)(N})"(N})" dx 

(50) 


The matrix functions are functions of a cross-sectional property which is shown as an argument in 
parentheses when appropriate, for example, Sfj(pA). The expressions for the nonzero terms in the mass and 
stiffness matrices given by equations (8) through (15) can now be expressed in terms of the matrix functions 
as follows: 


Kj = SfjipA) 

= Sfj(pA) + S} p (pl zz ) 

M?j = A(i)A(i) [slj(pA) + S/jVyy)] 
Mfj = Slj(plxx) 

= $J&A) 

KVj = S}J P (EI ZZ ) 


= X(i)XU)SlJ p (EI Y Y ) 
Kj = Sij(GJ) 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


where 



( i = 2 or 4) 
(Otherwise) 


(59) 


Explicit expressions for the five sets of matrix functions for t = l,oo and j = i,oc are given in appendix A. 
Expressions for j = 1, * — 1 are not given because the matrix functions are symmetric (i.e., S hJ = S Jt ). 


Element Mass and Stiffness Matrices 

The nonzero mass and stiffness terms given by equations (51) through (54) and equations (55) through 
(58), respectively, must be appropriately assembled to form the element mass and stiffness matrices. This 
procedure depends on the arrangement of the discrete degrees of freedom in the element vector of unknowns. 
The first 12 degrees of freedom in the vector of unknowns are the external degrees of freedom associated with 
the classical beam element. The higher order (internal) degrees of freedom are positioned after the external 
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degrees of freedom in the vector of unknowns. Specifically, all the internal q u, s are grouped together, then all 
the internal q v ' s are grouped together, and so on. The final arrangement of unknowns is given as follows: 


r <71 


r u(0) ' 

9 2 


u(0) 

<73 


u>(0) 

<7 4 


0(0) 

9 5 


— u/(0) 

<76 


</(0) 

97 


«(0 

98 


v(l) 

99 


w(l) 

910 


0(1) 

911 


-w'(l) 

912 


v'(l) 

913 


93 


1 < 


9l2+P„— 2 


ok 

9l3+P u — 2 



9l2+P M — 2+P„— 4 


ok 

9l3+P„— 2+P„— 4 



912+P u - 2+P„-4+P w -4 


ok 

9l3+P u 2+ P„ 4+Pm 4 


0 3 

k 9l2+P u -2+P„-4+P„,-4+P 0 -2 . 


> OPg . 


% 

I 

I 


Mass Matrix 

The element mass matrix is partitioned into several submatrices consistent with the arrangement of the 
degrees of freedom in the vector of unknowns given in equation (60). This leads to a matrix of the form 


M ee M Ij 

M|j 


1 

2 


M// 

0 

0 

0 



M h 

0 

0 

(61) 

sym 


Mfi 

0 





S 



The submatrix M ee is termed the external submatrix because it is associated only with the external degrees 
of freedom. This submatrix, which is given by the following equation, is a 12 by 12 matrix and reduces to the 
classical consistent mass matrix for constant cross-sectional properties: 
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M ee = 


Sf'l(flA) 0 


+Sj , ,\pJzz) 


0 tf 2 (pA) 0 

0 


S\ 2 (pA) 

+S\^(pl zz ) 


S}. s (M) 

+ S l^{p 1 zz) 


S\ a (pA) 

+Stf(ptrr) 


-SlJpA) 

-Stflp/yy) 


s%(pixx) 


syrn 


S} i2 (pA) 

+S&(plyy) 


0 


S}j,(pA) 


SjjpA) 

+Stf(pfzz) I 


+S]r(plyy) 

0 0 S^ 2 {plxx) 


-SljpA) 

-S\^(plyy) 


+s&Wi*) 


o 

S^(pA) 

+S\\ s {pi zz ) 


-S^JpA) 

-S£(plyy) 




«L|(M> 

+S£(pl zz ) 


0 

0 

Sl,(pA) 

+S&(plyy) 


+S‘>hv) 


S},i(pA) 

+S\'\{pl zz ) I 


s}. 4 (M) 

+S\1(plzz) I 


„ - 5 L(M) 

-SgUdyy) 

&2#(.PIxx) 0 

s| 4 (M) 

+S\'’(plyy) 


S\ JpA) 
*S\ p A (pl ZZ )\ 


(62) 

The submatrices and couple the external and internal degrees of freedom in the 

element mass matrix. These submatrices, which are shown in equations (63) through (66), are variable sized 
arrays of an order equal to 12 by the number of internal degrees of freedom corresponding to the particular 
continuous variable. For example, let P u represent the number of discrete degrees of freedom associated with 
u. Then, because there are always two external discrete degrees of freedom associated with u, the dimensions 
of M are 12 by P u — 2. It is of interest to note that the coupling terms associated with the higher order 
internal degrees of freedom become zero if a sufficient number of internal degrees of freedom are included. 
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S?, 4 (M) 

^1,5 (^) 

S?, 6 (M) 

o 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

o ... 

0 

0 

0 

0 

o ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

^2,3^^) 

S§ A (pA) 

*-*2,5 {pA) 

S^ipA) 

0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

o ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

o ... 

0 

0 

0 

0 

o .... 


0 

0 

0 

0 ..." 

s} i5 (M) + s\ p h (pi zz ) 

S\ 6 (pA) + S{P(pI zz ) ... 

5} 1 0 (M) + S\ p w (pl zz ) 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

S\z(pA) + s\r(pl zz ) 

S 2fi(pA) + S l ifiWzz) 

+ ^2 JO (P^ZZ) 

0 ... 

0 

0 

0 

0 ... 

S\ h {pA) + S\ p (pl zz ) 

Sls(pA) + SW(pI zz ) ... 

+ *^3J0 (plzz) 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

S\ b (pA) + S\ P b {pI zz ) 

S\^pA) + S\P{pI zz ) ... 

S\ M {pA) + S\%(pI zz ) 

0 ... 
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0 

0 

0 

0 ..." 

0 

0 

0 

0 ... 

5i\ 5 (pA) + 5}J(p/yy) 

+ ^i,6 (^y) 

5 wo(^) + s Mo(^y) 

0 ... 

0 

0 

0 

0 ... 

-5^(p/D-5g(p/ rK ) 

-S\^pA)-S\P{pI YY ) ... 

“*^2,10^^ “ ^2J0 (p?YY) 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

Sl 5 (pA) + Stf( p i yY ) 

S ls (pA) + s£(pWy) ... 

*^3,10 + *^3 JO (ptyy ) 

0 ... 

o 

0 ... 

0 

0 ... 

’co P 

t 

f 

-5} 6 (M) - Sl P 6 (pI YY ) ... 

-S\^(p a )~S\%(pI Y y) 

0 ... 

0 

0 

0 

0 ... 


0 

0 

0 

0 ..." 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

S l^ ! xx) 

S° lA (pixx) 

••• S u^ pl xx) 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 ... 

^Wxx) 

slMxx) 

••• S 2^p j xx) 

0 ... 

0 

0 

0 

0 ... 

0 

0 

0 

0 .... 


The mass submatrices corresponding to the internal degrees of freedom are given by the following equations 
(eqs. (67) through (70)). Note that these submatrices are symmetric and banded. Their size depends on the 
number of internal degrees of freedom employed for a particular continuous variable. Given that P v is the 
number of discrete degrees of freedom associated with v and that there are always four external degrees of 
freedom associated with v, the size of Mjj is P v — 4 by P v — 4. 


S? M (pA) Sfj 4 (pA ) 


$3.7 ^A) 0 0 


^(M) 


$4 8 (pA) 0 0 


M // = 


sym 


$ 5.5 b'4) ■‘^l M (pA) 


$$,<J (pA) 0 


(67) 
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S} r {pA) 

+4iwzz) 


M 


V _ 

11 - 


sym 


SL(PA) 

+Sl P (pIzz) 


Sin(pA) 

+Sl P n (pI zz ) 

0 

0 

SL(PA) 

s> 7 (m) 


Sl n {pA) 

0 

+Sl P (Phz) 

+S b P (P ! Zz) 


+ s l P uWzz) 


S)j(pA) 

Sl„(pA) 


Sfi 3 (M) 


+Sl l j(pI Z z) 

+S) p (plzz) 


+ S 1A’&(P I ZZ) 


(68) 


SljpA) 

+Sl1(plyy) 


M 


U) _ 
11 - 


sym 


Sl^pA) 

+Sl P (pIy Y ) 


Sin(pA) 

+Sl P u (pIyy) 

0 

0 

SU(PA) 

Sl 7 (pA) 


So.12<M) 

0 

+Sl p b (plyy) 

+Sl P (pIyy) 


+Sl P 2 (pIyy) 


S' 7J (pA) 

SljpA) 


S‘i 3 (^) 


+S' 7 p 7 lplyy) 

+S' 7 p (plyy) 


+S' 7 r u (plyy) 


( 69 ) 


M 


8 _ 
II ~ 


S 3,3(P I XX') S fl 4 (pI X x) ■■■ S^g(plxx) 
s° iA ( P ixx) S° 45 (pI XX ) 
sym ^(pfxx) S tfi( pl xx) 


0 0 0 

^S.loO^XX) 0 0 

ShMxx) o 


( 70 ) 


Stiffness Matrix 


The stiffness matrix is partitioned similar to the mass matrix and has the form 


K = 


K EE K'’ei K K», 

0 0 0 


K h 


K // 


0 

KW 


K 1 


II 


( 71 ) 
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The external stiffness submatrix T^ee is given by the 12 by 12 matrix in the following equation and reduces 
to the classical beam element stiffness matrix for constant cross-sectional properties. 

k ee = 


S^(EA) 0 0 0 

Sl P l p (EI zz ) 0 0 

S\ pp (EIyy) 0 -S\ pp (EIyy) 0 

S° P (GJ) 0 0 

S\ PP (EIyy) 0 


0 0 S^ P (EA) 0 0 

0 S\ PP (EI ZZ ) 0 S\ PP (EI ZZ ) 0 


0 

0 


0 0 
0 S\ PP (EI ZZ ) 


0 0 -S\ PP (EI YY ) 0 

0 0 S°£(GJ) 0 0 

0 -S l 2 PP (EIyy) 0 Sl PP (EIyy) 0 

0 0 S\ PP (EI ZZ ) 


sym 


S\ PP (EI ZZ ) 0 S\ PP (EI ZZ ) 0 

S% P (EA) 0 0 

S\ P /(Ehz) 0 

S\ pp (EIyy) 0 -S^ftEIyy) 0 

S% p 2 (GJ) 0 0 


0 0 
0 S\ PP (EI ZZ ) 


S\ PP (EIyy) 


S\ PP (EI ZZ ) 


(72) 


The submatrices K£j. K V EI , K and couple the external and internal degrees of freedom in the element 
stiffness matrix and are given by the following equations. These submatrices have the same dimensions and 
properties as the corresponding external-internal coupling mass submatrices (eqs. (63) through (66)) discussed 
previously. 


S^(EA) ea ) 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


\c u — 


0 0 0 

S^(EA) S% p 4 (EA) 0 

0 0 0 

0 0 0 

0 oo 

0 0 0 

0 0 0 


(73) 
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*EI 


TS-W _ 

*^EI ~ 




0 

0 

0 0 

" 



S\ P 5 P (EI ZZ ) 

S\ PP (EIzz) 

••• S\ PP (EI ZZ ) 0 





0 

0 

0 0 





0 

0 

0 0 





0 

0 

0 0 




S \ P f> P ( EI Zz) 

s \7^ EI zz) 

• s \7 {EI zz) o 



= 


0 

0 

0 0 




S$ P (EI ZZ ) 

$3,6 P (EIzz) 

• S\ P /(EI ZZ ) 0 





0 

0 

0 0 





0 

0 

0 0 





0 

0 

0 0 




S\ PP (EI ZZ ) 

S\ P /(EI ZZ ) 

... S\ PP {EI ZZ ) 0 



r 

0 


0 

0 

0 . 

1 


0 


0 

0 

0 ... 

S\ P 5 P (EI yy) 

5 ir( £/ yy) 

... 5}7(£;/ yr ) 

0 ... 


0 

0 

0 

0 ... 

- 

S\ PP (EI YY ) - 

-s£ p (ei yy ) 

... -S\ PP (E1 YY ) 

0 ... 


0 


0 

0 

o ... 


0 


0 

0 

o ... 


0 

0 

0 

0 ... 

s£ p (Ei YY ) 

Sl PP (EI YY ) 

... S\ PP {EI YY ) 

0 ... 


0 


0 

0 

0 ... 

- 

S\ PP (EI YY ) - 

-S\ PP (EI YY ) 

... -S\ PP (EI YY ) 

0 ... 

- 

0 

0 

0 

0 . 

• J 



0 

0 

0 0 ... - 




0 

0 

0 0 ... 




0 

0 

0 0 ... 




S°^(GJ) 

I S° X P (GJ) . 

. S» P (GJ) 0 ... 




0 

0 

0 0 ... 




0 

0 

0 0 ... 


K £/ - 







0 

0 

0 0 ... 




0 

0 

0 0 ... 




0 

0 

0 0 ... 




S% p ( GJ 

1 S$ P (GJ) . 

. S§£(GJ) 0 ... 




0 

0 

0 0 ... 




0 

0 

0 0 ... . 



( 74 ) 


( 75 ) 


( 76 ) 
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The stiffness submatrices associated with the internal degrees of freedom have similar characteristics as 
their mass counterparts and are defined as follows: 

S°$(FA) Sjfi(EA) Sgg(EA) 0 0 0 ... 

Sff(EA) S$ P (EA) S^ P (EA) 0 0 ... 

K u = (77) 

U sym Si” \(EA ) S°$(EA) S° P (EA) 0 ... 


Stf/( EI ZZ) S\ PP (EI ZZ ) ... S} PP (EI ZZ ) 0 n 0 

S^ P (EI ZZ ) S} PP (EI ZZ ) ... S\ PP (EI ZZ ) 0 0 

Stf P (EI ZZ ) S\ PP (EI ZZ ) ... S\ PP (EI ZZ ) 0 


S$ P Wyy) s& p (ei yy ) ■■■ Sl PP (EI YY ) 0 0 0 

S$ P (EI YY ) S$ p (EIyy) ••• S l 6 PP (EI YY ) 0 0 

sym S) pp (EI YY ) S\ pp (EI yy ) ■■■ S\ PP (EI YY ) 0 


&l P (GJ) S° P (GJ) ... S$ P (GJ) 0 0 0 

S° 4 P 4 (GJ) S° 4 p (GJ) ... S° 4 p (GJ) 0 0 

sym S% P (GJ) S$ P (GJ) ... S% P (GJ) 0 


All the relations needed to completely define the tapered p-version beam element for an arbitrary number 
of internal degrees of freedom are now in hand. These include the matrix definitions given in equations (61) 
through (80), the definition of the cross-sectional property polynomials given by equations (44) and (45), and 
the explicit expressions for the five sets of matrix functions given in appendix A. 

Numerical Validation and Preliminary Performance Analysis 

The beam element developed herein is capable of emulating four different types of beam elements: uniform 
h-version, uniform p-version, tapered h-version, and tapered p-version. A uniform (h- or p-version) element 
is created by restricting ai through 04 to be zero in the generic cross-sectional property polynomials. An 
h-version (uniform or tapered) element is created by restricting P u = Pg = 2 and P v = P w = 4. Once the 
uniform versions of the element are validated, they can be used to approximate a tapered geometry for the 
purpose of validating the tapered versions of the element. Similarly, the inversions of the element should 
converge to the same results as the h-versions of the element. 
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Computational Approach 

Numerical results for most of the cases analyzed were computed on an IBM Model 80 computer operating 
at 16 MHz and equipped with an 80387 math coprocessor and running PC-DOS 3.3. Because of computer 
core limitations, those uniform h-version cases containing 15 or more elements were run on a DEC 3200 
VAX workstation running VMS. However, all computations used eight-byte floating-point representation of 
real numbers conforming to the IEEE standard (default on the IBM computer but requires a /G compiler 
option on the VAX computer). Because the computer code used to generate these results was not written to 
take advantage of the symmetry and bandedness of the global mass and stiffness matrices, a fair comparison 
of compute times cannot be given. The only measure of computational efficiency presented is the number of 
degrees of freedom versus the accuracy of the models. The eigenvalue extraction technique used is a combination 
of the determinant search method and inverse iteration (ref. 16). The determinant search phase employs both 
the bisection method and the secant method to converge to an eigenvalue. Inverse iteration is then used to 
obtain the associated eigenvector. 


Validation of Uniform h- and p-Versions 

A uniform beam element can be validated rather easily because closed-form solutions to the axial and 
torsional vibrations and “nearly closed-form” solutions to the bending vibration problem exist. (See ref. 12.) 
The derivation of the closed-form expressions for these frequencies is presented in appendix B because the 
derivations commonly found in the literature do not include the rotary inertia terms. 

The uniform, cantilevered steel beam with a circular cross section shown in figure 3 was used to validate 
the uniform versions of the element developed herein. Representative h- and p- version finite-element models of 
the beam are also shown in figure 3. The material properties of the beam used in the numerical studies were 
as follows: 

E = 30.0 x 10 6 lb/in 2 


G = 11.6 x 10 6 lb/in 2 


0.284 

386.4 


lb-sec 2 /in 4 


For validation of the uniform versions of the finite element developed herein, the lowest frequencies of this 
beam were calculated “exactly” using the closed-form solutions in appendix B and numerically using both the 
h- and p-version elements. The first four bending frequencies, the first torsional frequency, and the first axial 
frequency, from appendix B, are 

f v i = 27.4651 Hz 


/„ 2 = 171.6155 Hz 


/„ 3 = 478.2670 Hz 


f vA = 930.8451 Hz 


foi = 523.4749 Hz 
/„ i = 841.7997 Hz 

These six frequencies will be regarded as “exact” for the purpose of validating the uniform h- and p-versions 
of the beam element. Table I shows the frequencies obtained using the uniform h-version of the finite element 
as the number of elements used was varied from 1 to 14. The exact frequencies are given at the bottom of 
the table for reference. Although only the first bending frequency truly converged with 14 elements, the other 
frequencies are nearly converged (within 0.06 percent). 

Table II shows the frequencies obtained using one p-version element as the number of internal degrees of 
freedom is increased from 0 to 13 for each of the four continuous displacements. That is, in table II, 13 internal 
degrees of freedom refers to 13 internal degrees of freedom each for u, v , w , and 0, for an aggregate totai of 
52 internal degrees of freedom. Again, the exact frequencies are noted at the bottom of the table for reference. 
It is clear that all the listed frequencies have converged with 13 internal degrees of freedom, and that the 
converged values are essentially equal to the exact values. In figures 4 through 9 the information in tables I 
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and II is plotted in terms of the percent error versus the number of degrees of freedom. The plots clearly show 
that both the p- version element and the h-version element converge to the exact answer, and that the p- version 
element has a faster convergence rate than the h-version element. 

Validation of Tapered h- and p- Versions 

The tapered steel beam with circular cross section shown in figure 10 was used to validate the tapered 
versions of the finite element developed herein. Also shown in figure 10 are representative h- and p-version 
models. Because the results shown in tables I and II have validated the uniform versions of the beam element, 
the exact frequencies of the tapered beam can be determined with a uniform h-version model if a sufficiently 
large number of elements are employed. 

The frequencies obtained for the tapered beam using uniform h-version finite elements are presented in 
table III. It should be noted that the frequencies in table III approach the converged values from below 
(whereas in tables I and II, the frequencies converge from above). Also, the number of elements required to 
reach convergence for the tapered beam is over seven times that required for a comparable uniform beam. 
Both these behaviors can be explained by the fact that, as uniform elements are added, the geometry of the 
model is changing. The geometry approaches a smoothly tapered beam as the number of elements in the model 
approaches infinity. Here, the model is assumed to be converged with 100 elements (600 degrees of freedom). 

The frequencies predicted by using a tapered h-version finite-element model are shown in table IV. In this 
case the geometry of the problem is represented exactly using one element, and therefore the geometry of the 
model is not changed as the number of elements is increased. It should be noted that now convergence is from 
above and is achieved using approximately one seventh as many elements as in table III. 

Table V shows the frequencies which were obtained for the tapered beam using one tapered p-version 
element and varying the number of internal degrees of freedom from 0 to 19. With 19 internal degrees of 
freedom, 22nd order polynomials are used for lateral bending and 20th order polynomials for the axial and 
torsional displacements. The reason for showing the results of such a high order element is to demonstrate the 
numerical stability of the element. It is seen that the p-version results in table V converge very close to the 
values predicted by the h-version beams in tables III and IV but with the use of far fewer degrees of freedom. 

The data in tables III through V are graphically depicted in figures 11 through 16 in terms of percent error 
versus the number of degrees of freedom. These figures show the dramatic improvement in convergence which 
is realized by using tapered h- and p-version elements instead of uniform h-version elements. 

Conclusions 

The derivation, and validation, of a new, tapered, p-version beam element which both facilitates convergence 
checks and produces a better convergence rate than nontapered, h-version beam elements has been described. 
These two characteristics complement each other and, when combined, provide a powerful and versatile beam 
element which is easy to use. The shape functions on which the element is based were derived by using 
orthogonality relationships which produce element matrices that are extremely well-conditioned and of a form 
allowing explicit integration in the derivation of the element matrices. The latter feature eliminates the need 
for numerical quadrature; thus, roundoff error is reduced. The shape functions are hierarchical such that higher 
order element matrices can use the element matrices from previous lower order analyses. This simplifies the 
derivation, coding, and validation of the element. The present form of the beam element has been derived in a 
manner which allows for an infinite number of internal degrees of freedom. The beam element has been tested 
with up to 22nd order C \- type shape functions and up to 20th order C°-type shape functions with no evidence 
of ill conditioning or significant roundoff error. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
June 23, 1989 
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Appendix A 

Expressions for Matrix Functions 

Explicit expressions for the matrix functions are presented in this appendix. Each matrix function is a 
function of the element length l and the coefficients from the generic cross-sectional property polynomials, oo 
through 04 . 

The definition of the 5° set is 


= J‘ P 4 (x)N?N* dx 


and the explicit integrations are 


0 /(4a 4 / 4 ^03/^ + 14a2/ 2 + 35 d\l + 140ag) 

^ = 420 


-,q /(10o 4 / 4 + 14o 3 / 3 -+■ 21a 2 / 2 + 35 Q.\l + 70ao) 

5l ’ 2 = “ 420 


c0 _ 

* 2,2 ~ 


/(6O04/ 4 + 70a 3 / 3 + 84 o 2 / 2 + 105aiZ + 140a 0 ) 


(Al) 

(A2) 

(A3) 

(A4) 


= -W2i - 3 [(? - 5)(33? 4 - 685? 3 + 4990? 2 - 149107 + 15492)a 4 / 4 
i=3 ’ 8 + 2 (? - 8)(23? 4 - 416? 3 + 2697? 2 - 7419? + 7305)a 3 / 3 

+ 6 (i - 8)(? - 7)(8? 3 - 90?' 2 + 322? - 365 )a 2 / 2 
+ 4 (? - 8 )(? - 7 )(* - 6)(i - 4)(4? - 5)oiZ 
+ 35(?' - 8)(? - 7 )(i - 6)(i - 5)(2? - 5)a 0 ] 

/[168(587? 5 - 15825? 4 + 166115? 3 - 846075? 2 + 2087798? - 1994340)] 

S° hi =0 

7=9,00 

5^ = lV2i^3 [(283? 5 - 7825? 4 + 84335? 3 - 442025? 2 + 1124532? - 1108260)a 4 / 4 
i=3 ’ 8 + 5 (? - 8)(50? 4 - 995? 3 + 7224? 2 - 22641? -f 25770 )a 3 / 3 

+ 3 (? - 8)(? - 7)(68? 3 - 905? 2 4- 3907? - 5440)a 2 7 2 
+ 3(? - 8)(? - 7)(? - 6)(52? 2 - 413? + 785)oii 
+ 35(? - 8)(? - 7 )(? - 6)(? - 5)(4? - 13)a 0 ] 

/[168(587? 5 - 15825? 4 + 166115? 3 - 846075? 2 + 2087798? - 1994340)] 

sl =0 

7=9,00 

= 7(2? - 3) [3(7? 4 - 42? 3 + 7 ? 2 + 168? + 60)a 4 / 4 
*=3,oo + ( 2 j + 3 )(2? - 9)(7? 2 - 21? - 10 )o 3 7 3 

+ 2(2? + 3)(2? - 9)(5? 2 - 15? - 8 )a 2 / 2 
+ 4(2? + 3)(2? + 1)(2? - 7)(2 i - 9)(oi l + 2a 0 )] 

/[16(2? + 3)(2? + 1)(2? - 1)(2? - 3)(2? - 5)(2? - 7)(2? - 9)] 


(A5) 

(A6) 


(A7) 

(A8) 
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(A10) 


Sj i+1 = -l 2 VW^3V2T^l[6(i - 1)(* 2 -2 * - 5 )a 4 Z 3 

i=3 ’°° + (* - 1)(7* 2 - 14* - 36)a 3 Z 2 + 2(i - 1)(2* + 3)(2* - 7 )(a 2 / + a x )] 

/[16(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7)] 

s i,i+2 = -lV2i~^3\/Ti + 1 [15(* 2 - * - 10)(* 2 - * - 8)a 4 Z 4 

i='3,oo + 8 (j + 2 )(i - 3)(2* + 5)(2* - 7)a 3 / 3 + 8(2 i + 5)(2 * - 7)(2* 2 - 2* - 9)a 2 Z 2 

+ 8(2* 4- 5)(2 * + 3) (2* - 5) (2* - 7)(a x / + 2a 0 )] 

/[64(2* + 5)(2* + 3)(2 * + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7)] 

Sg i+3 = l 2 \/2i — 3\/2* + 3[l0*(* 2 - 7)a 4 / 3 + 7(13* 2 - 85)a 3 Z 2 
i=3, (X) +4*(2* + 5)(2* — 5)(a 2 Z + a x )] 

/[32(2t + 5)(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5)] 

Sj i+4 = -Z 3 v/2T^3v / 2iT5[(* + 1)*(13* 2 + 13/ - 116)a 4 Z 2 
i=3,oo + (j + l)i(2* + 7)(2 * - 5)(3 a 3 Z + 2a 2 )] 

/[32(2t + 7)(2 * + 5)(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5)] 

g / 4 \/2* — 3\/2* 4- 7(* + 2 )(* + l)*(2a 4 Z 4- a 3 ) 

= 32(2* + 7)(2* 4- 5)(2* + 3)(2t + 1)(2* - 1)(2* - 3) 

g —\fli — 3^2* + 9 (* 4- 3)(* + 2 )(* 4~ l)*a 4 / 3 

>* +6 = 64(2* 4- 9)(2* + 7)(2* 4- 5)(2* + 3)(2* 4- 1)(2* - 1)(2* - 3) 

cQ . . — Q 

°l,l+3 — u 

7—3,00, j=7,oo 

The definition, of the S op set is 

S?X = j‘p*(x)(N?)'(N"y dx 

and the explicit integrations are 

q p 12a 4 / 4 + 15a 3 / 3 + 20a 2 Z 2 + 30a x Z 4~ 60ao 

1.1 ~ 60/ 

q p 12 a 4 / 4 + 15a 3 / 3 + 20 a 2 / 2 4~ 30a x / + 60ao 

1.2 = 60/ 

gp 12a 4 / 4 + 15a 3 / 3 4~ 20a 2 / 2 4“ 30a x Z + 60a 3 

2,2 = 60/ 

SfJ = x/2r03[3(29* 3 - 406* 2 4- 1833* - 2644)a 4 / 3 

i=3 ’ 6 4- 3(* - 6)(25* 2 - 202* 4- 387)a 3 / 2 4- 2(* - 6)(* - 5)(26* - 83)a 2 / 

4- 10(t - 6)(* - 5)(* - 4)ai]/[30(47* 3 - 628* 2 4- 2729* - 3792)] 


(All) 

(M2) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 
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Sff = 0 

2 = 7,00 


(A22) 


= -y/2 i - 3[3(29* 3 - 406* 2 + 1833* - 2644 )a 4 / 3 + 3 (* - 6)(25* 2 - 202* + 387 )a 3 / 2 
i=3,e + 2 (i - 6 )(* - 5)(26* - 83)a 2 / + 10(t - 6)(i - 5 )(* - 4)ai] 

/ [30(47« 3 - 628* 2 + 2729? - 3792)] ( A23) 


=0 

2=7,00 

(A24) 

sff = [(35 i 4 - 210* 3 + 305* 2 + 30 i - 88) (2 * - 3 )a 4 / 4 
*=3,oo + 2(2 i + l)(2 * - 3)(2 i - 7)(5 i 2 - 15* + 7 )a 3 / 3 

+ 4(2* + 1)(2* - 3)(2* - 7)(3* 2 - 9* + 4)a 2 / 2 
+ 4(2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7){a x l + 2a 0 )] 
/ [8/(2* + 1)(2* - 1)(2* - 3)(2t - 5)(2* - 7)] 

(A25) 

S£f +1 = — \/2* — 3\/2* — 1 [2(* - 1)(7* 2 - 14* - 8)a 4 Z 3 
i=3,oo + 3 (j _ !)( 5 j 2 _ 1Q - _ 6 ) a3/ 2 

+ 4(t - 1)(2* + 1)(2* - 5)(o 2 / + ai)] 

/ [8(2* + 1)(2* - 1)(2* - 3)(2* - 5)] 

(A26) 

S i,f+2 = \/2* - 3>/2 * + 1[2*(* - 1)(7* 2 -7 * - 26)a 4 Z 3 
i=3 ’°° + *(* - 1)(2* + 3)(2* - 5)(3a 3 / 2 + 2a 2 /)] 

/ [8(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2t - 5)] 

(A27) 

op — 3\/2i -|- 3(i + 1 )i(i — l)(2a4/ 3 + <23 l 2 ) 

is £2 ~ 8(2* + 3)(2* + 1)(2* - 1)(2* - 3) 

(A28) 

cOP yj2i — 3\/2 i 5(z “I- 2)(i -b 1 )i{i — 1)g-4 / 3 

> 3 ^ “ 16(2* + 5)(2* + 3 ) (2* + 1)(2* - 1)(2* - 3) 

(A29) 

cQP — n 

^ 1 , 1 +j — u 

i=3,oo, j=5,oo 

(A30) 


The definition of the S 1 set is 

Slj = j‘p 2 (x)NlNj dx (A31) 


I 
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and the explicit integrations are 


Z(19a2^ 2 + 54a\l -b 234ao) 
6l > 4 “ 630 

(A32) 

ol l 2 {17a 2 l 2 + 42aiZ + 132a 0 ) 

1)2 “ 2520 

(A33) 

Z(46a 2 Z 2 + 81ajZ + 162ao) 
1)3 “ 1260 

(A34) 

Z 2 (19a 2 Z 2 + 36a iZ + 78ao) 
1,4 “ 2520 

(A35) 

Z 3 (4a 2 Z 2 + 9 a\l 4- 24ao) 
2 ’ 2 “ 2520 

(A36) 

Z 2 (25a 2 Z 2 + 42aiZ + 78ao) 
2 ’ 3 ~ 2520 

(A37) 

Z 3 (5 a 2 Z 2 + 9ajZ + 18ao) 
2 ’ 4 “ 2520 

(A38) 

el _ Z(145a 2 Z 2 + 180aiZ + 234a 0 ) 
3 ’ 3 “ 630 

(A39) 

c l Z 2 (65a 2 Z 2 4- 90ajZ 4- 132ao) 

3 ’ 4 “ 2520 

(A40) 

ol _ Z 3 (10a 2 Z 2 + 15aiZ + 24a 0 ) 
4 ’ 4 ~ 2520 

(A41) 


5j j = v/2*-5[(206* 5 - 7725 * 4 + 114020* 3 - 827175* 2 + 2947034* - 4122960)a 2 Z 3 
*= 5 ’ 10 + (i - 10)(173 *' 4 - 5005 * 3 + 53215 * 2 - 245915* + 416292)aiZ 2 
+ 2 (*' - 10) (i - 9 )(* - 7) (69 *' 2 - 1136* + 3892)a 0 Z] 

/[1260(4759* 5 - 174190* 4 + 2512445* 3 - 178441 10* 2 + 62404896* - 85980000)] (A42) 

Si = 0 (A43) 

7 = 11,00 

Si = \/2* — 5[(62* 5 - 2315* 4 + 34000* 3 - 245305* 2 + 868878* - 1208400)a 2 Z 4 
*=5,10 + 2 (j _ io)(* - 8)(* - 6)(33* 2 - 468* + 1499)ai/ 3 

+ (* - 10)(* - 9)(26* 3 - 66 1* 2 + 4961* - 11476)a 0 Z 2 ] 

/ [10080(434* 5 - 15875* 4 + 229060* 3 - 1628305* 2 + 5701086* - 7864620)] ( A44) 


Sl i = — y/zT—l) [5( 156* 5 - 5983* 4 + 90554* 3 - 675437 * 2 + 2480302* - 3583344)a 2 Z 3 
*=5,lO + _ io)(l09i 4 - 3167* 3 + 33983* 2 - 159337* + 274980)ai* 2 

+ 2 (* - 10)(* - 9)(* - 7) (1 ll* 2 - 1724* + 5908)a 0 /] 

/[l260(4759i 5 - 174190* 4 + 2512445* 3 - 17844110* 2 + 62404896* - 85980000)] 

Sli = 0 
1=11,00 

s\ ti = -V2f^5[Q{8i 5 - 275* 4 + 3670* 3 - 23685* 2 + 73642* - 87900)a 2 / 4 
*=5,10 + 5 (j _ 10)(22* 4 - 573* 3 + 5474* 2 - 22743* + 34716)oi/ 3 

+ (* - 10)(* - 9)(176* 3 - 3141* 2 + 18371* - 35276)a 0 / 2 ] 

/[10080(434* 5 - 15875* 4 + 229060* 3 - 1628305* 2 + 5701086* - 7864620)] 

sj,i =0 

7=11, OC 

x _ (2* - 5)/[(7* 2 - 35* - 18)a 2 / 2 + 3(2* + 1)(2* - ll)(ai* + 2 a 0 )] 
i-i'L ~ 16(2 * + 1)(2* - 1)(2* - 3)(2* - 5) (2* - 7)(2* - 9)(2* - 11) 


j \/2 1 — 5\/2 * — 3(* — 2 )* 2 (a 2 * + a\) 

>5^ = ” 16(2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7)(2* - 9) 

, _ *\/2*^5\/2z^T[( 17* 2 — 51* — 126)a 2 Z 2 4- 8(2* + 3)(2* — 9)(ai* + 2a 0 )] 

m+ 2 “ _ 64(2* + 3)(2* + 1) (2* - 1)(2* - 3)(2* - 5)(2* - 7)(2* - 9) 

l — 0,00 

x _ 3/ 2 \/2* - 5\/2 * + 1 (* - l)(a 2 / + ai) 

it ££ " 32(2* + 3) (2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7) 

j _ l\/2i — 5\/2 * + 3[(* + 4 )(* - 5 )a 2 * 2 + (2* + 5)(2* - 7)(a\l + 2ao)] 

> i+4 = 32(2* + 5)(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7) 


x _ / 2 \/2* - 5^2* + 5*(a 2 * + oi) 

= _ 32(2* + 5) (2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5) 

j \/2* — 5\/2* + 7 (* + 1 )*a 2 Z 3 

= 64(2* + 7)(2* + 5)(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2* - 5) 

S h +j = o 

5,oo, j=7,oo 


(A46) 

(A47) 

(A48) 

(A49) 

(A50) 

(A51) 

(A52) 

(A53) 

(A54) 

(A55) 

(A56) 

(A57) 
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The definition of the S lP set is 


S}P = J‘ P,(x)(Nl)'(Njy dx 
and the explicit integrations are 


plP 10(24 “h 15(23/ 3 4 24a 2 / 2 4 42 a\l 4* 84og 
s hl = 


70/ 


j p 5 C 24/ 4 *4" 4 10(22 4" 14(2^/ 4 14(2q 

140 


cir _ 
^1,2 “ 


ip IOG 4/ 4 T 16 ( 23 /^ T 24a 2 / 2 4 42ai/ 4 84(2o 

70/ 


oiP _ 
^1,3 ” 


IP _ 5a4/ 4 4 503/^ 4 4a 2 / 2 — 14ag 

~~ 140 


Qir _ _ 
^1,4 “ 


C 2 P /( 8 C 24/ 4 4 1 1 a. 3 /^ 4 16a 2 / 2 4 28(2 j/ 4 112ao) 
* 42,2 = 


840 


^P 5 « 4/ 4 4 7 ( 23 /^ T 10a 2 / 2 4- 14(2^/ 4- 14<2q 


^3-- 140 


'rlP /(IOG 4/ 4 4* II 03/ 3 4- 12a 2 / 2 4" 14oi/ 4" 28flo) 


cl/' _ ' 

^ 4 ~~ 840 


^ 2 P 1 0c2 4 Z 4 4 I 503 Z 3 4- 24a 2 Z 2 4* 42(2^/ 4" 84(2 q) 


cir _ 

3 ’ 3 70 / 


ol p __ 5a4 / 4 4 5a3 / 3 4 4a 2 / 2 — 14ag 
34 = 140 


IP /( 6 OG 4 Z 4 4 65 ( 23 / 3 4 72a 2 Z 2 4 84(2 1 Z 4 112ao) 

840 


*^ 4,4 “ 


5 }J = v / 2T 3 5[12(36/ 5 - 1330/ 4 4 19340/ 3 - 138305/ 2 4 4863192 - 672680)a 4 / 4 

i =$ A 0 + 1Q ^ _ 10 ^ _ 6 )( 47i 3 _ 993^2 + 68 q 3 - __ 15053 )a 3/ 3 

4 (i — 10)(t - 9)(419Z 3 - 7744Z 2 4 471292 - 94384)a 2 / 2 
4 (2 - 10)(i - 9)(i - 8 )( 269/ 2 - 3112* 4 8823)ai/ 

4 330(2 - 10) (2 - 9)(* - 8)(* - 7)(* - 5)a 0 ] 

/[840/(687z 5 - 25346Z 4 4 368441/ 3 - 2636686/ 2 4 9287944/ - 12883260)] 

stf = o 

i= ll,oo 


(A58) 

(A59) 

(A60) 

(A61) 

(A62) 

(A63) 

(A64) 

(A65) 

(A66) 

(A67) 

(A68) 


(A69) 

(A70) 
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= V2i - 5[8(103* 5 - 3795* 4 + 55020* 3 - 392190/ 2 + 1374332* - 1894305)a 4 / 4 
*= 5,10 + (j _ io)(991* 4 - 26735* 3 + 266005* 2 - 1157585* + 1860084)a 3 / 3 

+ 4 (* - 10)(* - 9)(i - 6)(247* 2 - 3010* + 8857 )a 2 / 2 
+ 3 (* - 10)(* - 9)(* - 8)(159* 2 - 1792* + 4973)<*i Z 
+ 6 (*' - 10)(* - 9)(* - 8)(* - 7)(158* - 783)a 0 ] 

/[1680(1918* 5 - 71033* 4 + 1035968* 3 - 7434583* 2 + 26251650* - 36487620)] 

stf = o 
2=11,00 

= — \/2* — 5[l2(36* 5 - 1330* 4 + 19340i 3 - 138305/ 2 + 486319* - 672680)a 4 / 4 
i=5, 10 + 10 (i _ io)(i - 6)(47* 3 - 993* 2 + 6803* - 15053)a 3 / 3 

+ (* - 10)(* - 9)(419* 3 - 7744* 2 + 47129* - 94384)a 2 / 2 
+ (* - 10)(* - 9)(* - 8)(269* 2 - 3112* + 8823)ai/ 

+ 330(* - 10)(* - 9)(* - 8)(Z - 7)(* - 5)a 0 ] 

/[840/(687* 5 - 25346* 4 + 368441* 3 - 2636686* 2 + 9287944/ - 12883260)] 

-o 

2=11,00 

= \/2*-5[2(1412* 5 - 52945* 4 + 782650* 3 - 5698355Z 2 + 20426118/ - 28825620)a 4 Z 4 
*=5,10 + 1Q ^ _ 10 )(247* 4 - 6828* 3 + 69745/ 2 - 311950/ + 515286)a 3 / 3 

+ 2 (/ - 10)(* - 9)(949/ 3 - 17989* 2 + 112234* - 230124)a 2 / 2 
+ 3(* - 10)(* - 9)(i - 8)(393* 2 - 4614* + 13261)ai/ 

+ 6(Z - 10)(Z - 9)(Z - 8)(* - 7)(172Z - 867)a 0 ] 

/[1680(1918* 5 - 71033* 4 + 1035968* 3 - 7434583Z 2 + 26251650/ - 36487620)] 

C =o 

2 = 11,00 

Sjf = (2 i - 5)[3(7* 4 - 70Z 3 + 175Z 2 - 52)a 4 Z 4 + (2 i + 1)(2 * - 11)(7Z 2 - 35* + 18)a 3 / 3 
*=5,oo + 2 ( 2 i + 1)(2* - 11)(5Z 2 - 25* + 12)o 2 Z 2 

+ 4(2* + 1)(2* - 1)(2* - 9)(2Z - ll)(oiZ + 2a 0 )] 

/ [16/(2* + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7)(2 * - 9)(2* - 11)] 

Sjf +1 = -y/2i^5^/2i^3[6{i - 2)(Z 2 - 4* - 2)a 4 Z 3 + (Z - 2)(7* 2 - 28* - 15)a 3 Z 2 
*=5,oo + 2 (/ - 2)(2 * + 1)(2* - 9)(a 2 / + «i)] 

/[16(2* + 1)(2Z - 1)(2* - 3)(2* - 5)(2Z - 7)(2Z - 9)] 

S}f + 2 = - 3* - 8)(/ 2 - 3/ - 6)a 4 / 4 + 8 (* + 1)(* - 4)(2Z + 3)(2 i - 9 )a 3 / 3 

*=5,°c + 8 ( 2 - + 3 )( 2 - _ 9 )( 2i 2 _ 6i _ 5 ) ffl2/ 2 + 8 ( 2i + 3)( 2 j + 1)(2* - 7)(2 * - 9 )(a 4 / + 2 a 0 )] 

/ [64/(2* + 3)(2Z + 1)(2Z - 1)(2* - 3)(2* - 5)(2Z - 7)(2* - 9)] 


(A71) 

(A72) 


(A73) 

(A74) 


(A75) 

(A76) 

(A77) 

(A78) 

(A79) 
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(A80) 


Si,i + 3 = v/2T^/2T+T[lO(i - 1)(* 2 - 2 i - 6 )a 4 Z 3 + (i - l)(13i 2 - 26* - 72)a 3 Z 2 

i=5,oo + 4 (j _ i)( 2 » + 3 )( 2 i - 7 )(a 2 / + aj)] 

/ [32(2* + 3)(2i + 1)(2* - 1)(2* - 3)(2* - 5)(2* - 7)] 

S} p + 4 = -v / 2i-5v / 2i + 3[i(i - 1)(13* 2 - 13 i - 116 )a 4 / 3 
i=5 ’°° + - 1) (2* + 5)(2* - 7)(3a 3 Z 2 + 2 a 2 Z)] 

/ [32(2* + 5)(2* + 3)(2i + l)(2i - 1)(2* - 3)(2t - 5)(2t - 7)] 

p\P \/2i — 5 v^2* + 5 (i + l)i(i — l)(2a4/ 3 + a 3 Z 2 ) 

= 32(2* + 5)(2* + 3)(2* + 1)(2Z - 1)(2Z - 3)(2i - 5) 

p x/2* - 5\/2t 4- 7(? 4- 2 )(j + 1)x(j 1 ) 1 X 4/ 3 

*5 ? +6 = — 64(2j + 7)(2» + 5)(2* + 3)(2* + 1)(2* - 1)(2* - 3)(2t - 5) 


cl P . 

s t,t+j 

7=5, oc, j=7,oc 


= 0 


Finally, the definition of the S 1PP set is 

S\ pp = ^P a {x){N})"{N})" dx 

and the explicit integrations are 


\pp 3 ( 4404 /^ + 49a 3 Z 3 + 56a2/ 2 + 70oiZ + 140ao) 


c nr _ 

* 1,1 — 


357 3 


ipp 3804 /^ + 42 o 3 / 3 + 49d2 / 2 + 70ai/ + 210oo 


oirr _ 

* 1,2 - 


35Z 2 


5 


ipp _ 
1,3 - 


3(44a 4 Z 4 + 49a 3 Z 3 + 56a 2 Z 2 + 70aif + 140o 0 ) 
35Z 3 


clPP - 

*1,4 — 


clPP - 

* 2,2 ~ 


clPP _ 
*2,3 - 


clPP _ 
*2,4 - 


ClPP - 
*3,3 “ 


94a4Z 4 + 105a 3 Z 3 + 119a2/ 2 + 140aiZ + 210ao 
35 Z 2 

36a 4 Z 4 + 42a 3 Z 3 + 56a 2 Z 2 + 105aiZ + 420a 0 
105Z 

38a 4 Z 4 + 42a 3 Z 3 + 49a 2 Z 2 + 70aiZ + 210a 0 
35Z 2 

78 a 4 Z 4 + 84a 3 Z 3 + 91a 2 Z 2 + 105aiZ + 210a 0 
105Z 

3(44a 4 Z 4 + 49a 3 Z 3 + 56a 2 Z 2 + 70a x l + 140a 0 ) 
35Z 3 


(A81) 

(A82) 

(A83) 

(A84) 

(A85) 

(A86) 

(A87) 

(A88) 

(A89) 

(A90) 

(A91) 

(A92) 

(A93) 
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olP P 94 <z 4 / 4 + 105 a 3 / 3 + 119a 2 / 2 + 140ai/ + 210a 0 

53 ’ 4 35/2 

(A94) 

olPP 204a 4 / 4 + 231a 3 / 3 + 266a 2 / 2 + 315ai/ + 420a 0 
^ 4 ’ 4 105/ 

(A95) 

S\$ p = -\/2t-5[3(45 Z 3 - 904Z 2 + 5965/ - 12902)a 4 / 3 
i=5 ’ 8 + 3(Z - 8)(38Z 2 - 463/ + 1383)a 3 / 2 + 2(Z - 8 )(/ - 7)(41Z - 219)a 2 / 

+ 28(i - 8 )(/ - 7 )(/ - 6 )oi]/[7/ 2 (68/ 3 - 1329Z 2 + 8551/ - 17970)] 

(A96) 

o 

II 

oT.ii 

(A97) 

S\ P { P = V2i - 5[2(50/ 3 - 996Z 2 + 6526/ - 14061 )a 4 / 3 
i=5,S + 12 (i _ 8)(7* 2 - 85/ + 255)a 3 / 2 + (/ - 8)(/ - 7)(62Z - 345)a 2 Z 

+ 42(i - 8)(Z - 7)(i - 6)oi]/[252Z(4/ 3 - 77 / 2 + 483/ - 995)] 

(A98) 

o 

i=9,oo 

(A99) 

S$ P = v/2/-5[3(45/ 3 - 904/ 2 + 5965/ - 12902)a 4 / 3 
i=5 ’ 8 + 3(i - 8)(38/ 2 - 463/ + 1383)a 3 / 2 + 2 (i - 8)(Z - 7)(41Z - 219)o 2 / 

+ 28(/ - 8)(/ - 7 )(/ - 6)ai]/[7/ 2 (68/ 3 - 1329/ 2 + 8551/ - 17970)] 

(A100) 

8 

II 

o 

(A101) 

S\? P = V2/-5[2(74/ 3 - 1485/ 2 + 9811/ - 21333)a 4 / 3 
i=5 ’ 8 + 3(Z - 8) (2/ - ll)(19i - 129)a 3 / 2 + (/ - 8)(/ - 7)(76/ - 429)a 2 Z 

+ 42(i - 8)(/ - 7 )(/ - 6)oi]/[252/(4/ 3 - 77 Z 2 + 483/ - 995)] 

(A102) 

sir = 0 

i= 9,oo 

(A103) 

slf p = ( 2 / - 5) [(35/ 4 - 350/ 3 + 1145Z 2 - 1350/ + 432)a 4 / 4 

i=5,oo + 2 ^ 2i _ 2i _ 9 )( 5i 2 _ 25i + 27)a 3 / 3 + 4(2/ - 1)(2Z - 9)(3/ 2 - 15/ + 16)a 2 / 2 
+ 4(2/ - 1)(2Z - 3) (2/ - 7) (2/ - 9)(a 4 Z + 2 a 0 )] 

/ [8/ 3 (2Z - 1)(2Z - 3)(2Z - 5)(2Z - 7)(2 / - 9)] 

(A104) 

S li+l = -\/2 / — 5^/2/ — 3[2(Z - 2)(7Z 2 - 28/ + 13)a 4 / 3 + 3(Z - 2)(5/ 2 - 20/ + 9)a 3 / 2 

i= 5 ,oo + _ 2)(2i - 1)(2 / - 7 )(a 2 Z + ai )]/[8Z 2 (2Z - 1)(2Z - 3)(2Z - 5)(2Z - 7)] 

(A105) 

S li+2 = V2Z - 5^2/ - 1 [2(Z - 1)(Z - 2)(7/ 2 - 21/ - 12)a 4 / 2 
i=5,oo + (j_ !)(< _ 2) (2/ + 1)(2Z - 7)(3a 3 / + 2a 2 )] 



/[8/(2Z + 1)(2 Z - 1)(2Z - 3)(2 i - 5)(2 i - 7)] (A106) 
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(A107) 

(A108) 

(A109) 
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Appendix B 

Derivation of Vibrations of Uniform Beam 

The derivation of the closed-form solutions to the vibrations of a uniform beam is presented in this appendix. 
First the bending vibration problem is solved, and then the torsional and axial vibration problems are solved. 


Consider the translational displacement v only. Using the expressions for the kinetic and strain energies 
given in equations (2) and (3), respectively, and assuming harmonic motion give the following governing 
differential equation: 


v m +e&v" p A‘ , *v = o 

h EIzz 

When taking into account the cantilever conditions, the geometric boundary conditions are 

(Bl) 


<2 

II 

O 

(B2) 


O 

II 

O 

(B3) 

and the natural boundary conditions are 

v"(l) = 0 

(B4) 

»"'(() + ^»'(0 = 0 

(B5) 

The general solution to equation (Bl) is 



v{x) = ki e blX + k 2 

e blX + k% cos(b 2 x) + £4 sin(&2x) 

(B6) 

where k\, k 2 , £3, and k± are constants that depend upon the boundary conditions, and b\ and b 2 

by 

are defined 

= \ 

tJa^ + aAo - a 2 

(B7) 

2 

h = \ 

^JaI + aAq + a 2 
2 

(B8) 

where 

, p a ^I 

° EI ZZ 

(B9) 


to 

II 

H3L 

(BIO) 


Substituting the general solution given by equation (B6) into the boundary conditions (eqs. (B2) through (B5)) 
results in four linear equations for k\ through k 4 that can be expressed in matrix form as follows: 


1 1 1 0 1 A* ' 

b\ —b\ 0 b 2 k 2 

b\ e bl1 b\ e ~ bl1 -&2Cos(M) < 

lbi(bl + A 2 ) e bl1 -bi(b'( + A 2 ) e ~° 1 l1 b 2 (bl ~ ^2)sin(6 2 0 -02(^2 - A 2 )cos(u2l)j \ -4 , 



The bending natural frequencies of the beam are those which make the determinant of the 4 by 4 matrix 
in equation (Bll) equal to zero. For the example described in the main text of this report, the first four 
frequencies are 
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fvl = 27.4651 Hz 
f„2 = 171.6155 Hz 
f v 3 = 478.2670 Hz 
f V 4 = 930.8451 Hz 

Now consider only the torsional displacement 6. Taking equations (2) and (3) as the kinetic and strain 
energies, respectively, and assuming harmonic motion give the following governing differential equation: 



s" + l,lx /f»e = o 

GJ 

(B12) 

The boundary conditions are 

a/fO — n 

is y \J ) — u 

(B13) 


e'(i) = o 

(B14) 

With b 2 = PlX Qy ^ , the general solution to equation (B12) can be written as 



6(x) = k\ cos(fcr) + &2 sin (bx) 

(B15) 


Substituting equation (B15) into the boundary conditions leads to the 2 by 2 matrix system of linear equations: 

1 0 ' 

_ — sin(6/) cos(6f) . 

Setting the determinant equal to zero and solving for ujq yields 

(n = 1,2, . . .) (B17) 

For the beam of figure 3, the first torsional frequency is 

f ei = 523.4749 Hz 

The governing differential equation for the axial displacement of the beam is of the same form as that for 
the torsional displacement given in equation (B12). Only the constant b is different. The natural frequencies 
can be shown to be 

(n = 1, 2, . . .) (B18) 

For the beam of figure 3, the first axial frequency is 

f u 1 = 841.7997 Hz 



nn I ~GJ 

Ue = 1 2l\J^ 


£-! 


(B16) 


30 


References 


1. Cronkhite, J. D.; Berry, V. L.; and Dompka, R. V.: Summary of the Modeling and Test Correlations of a NASTRAN Finite 
Element Vibrations Model for the AH-1G Helicopter. NASA CR- 178201, 1987. 

2. Gabel, R.; Reed, D.; Ricks, R.; and Kesack, W.: Planning, Creating and Documenting a NASTRAN Finite Element Model 
of a Modern Helicopter. Rotorcraft Dynamics 1984, NASA CP-2400, 1985, pp. 307 323. 

3. Toossi, Mostafa; Straub, Friedrich; Weisenburger, Richard; and Ferg, Douglas: AH-64A Shake Test and Correlation With 
NASTRAN Finite Element Model Predictions. 43rd Annual Forum Proceedings, Volume I, American Helicopter Soc., c.1987, 
pp. 197-204. 

4. Durno, Jason A.; Howland, Guy R.; and Twomey, William J.: Comparison of Black Hawk Shake Test Results With NASTRAN 
Finite Element Analysis. 43rd Annual Forum Proceedings, Volume II, American Helicopter Soc., c.1987, pp. 511-526. 

5. Lindberg, G. M.: Vibration of Non-Uniform Beams. Aeronaut. Q., vol. XIV, pt. 4, Nov. 1963, pp. 387-395. 

6. Sanger, D. J.: Transverse Vibration of a Class of Non-Uniform Beams. J. Mech. Eng . Sci., vol. 10, no. 2, Apr. 1968, 

pp. 111-120. 

7. Thomas, J.; and Dokumaci, E.: Improved Finite Elements for Vibration Analysis of Tapered Beams. Aeronaut. Q., vol. 24, 
pt. 1, Feb. 1973, pp. 39-46. 

8. To, C. W. S.: Higher Order Tapered Beam Finite Elements for Vibration Analysis. J. Sound & Vibration, vol. 63, no. 1, 
Mar. 8, 1979, pp. 33-50. 

9. To, C. W. S.: A Linearly Tapered Beam Finite Element Incorporating Shear Deformation and Rotary Inertia for Vibration 
Analysis. J. Sound & Vibration, vol. 78, no. 4, Oct. 22, 1981, pp. 475-484. 

10. Hodges, Dewey H.; Hopkins, A. Stewart; Kunz, Donald L.; and Hinnant, Howard E.: Introduction to GRASP — General 
Rotorcraft Aeromechanical Stability Program — A Modern Approach to Rotorcraft Modeling. J. American Helicopter Soc., 
vol. 32, no. 2, Apr. 1987, pp. 78-90. 

11. Craig, Roy R., Jr.: Structural Dynamics — An Introduction to Computer Methods. John Wiley & Sons, Inc., c.1981, 
pp. 198-199. 

12. Hurty, Walter C.; and Rubinstein, Moshe F.: Dynamics of Structures. Prentice-Hall, Inc., c.1964, pp. 177-201. 

13. Rivello, Robert M.: Theory and Analysis of Flight Structures. McGraw-Hill, Inc., c.1969, pp. 168 and 215. 

14. MSC/NASTRAN User's Manual— MSC/ NASTRAN Version 65. MSR-39, MacNeal-Schwendler Corp., Nov. 1985. 

15. Gallagher, Richard H.: Finite Element Analysis— Fundamentals. Prentice- Hall, Inc., c.1975, pp. Ill, 115, and 135-143. 

16. Bathe, Klaus- Jurgen: Finite Element Procedures in Engineering Analysis. Prentice- Hall, Inc., c.1982, pp. 557-672. 


31 



Table I. Convergence of Frequencies for Uniform, h- Version Beam 


fv \i Hz 

fv 2, Hz 

fv 3. Hz 

fv 4i Hz 

fe l, Hz 

fu l, Hz 

27.5953 

27.4784 

270.3419 

173.0589 

581.7178 

1649.2712 

577.2137 

537.0115 

928.2170 

863.5679 

27.4679 

172.1745 

484.1519 

1081.3530 

529.4733 

851.4457 

27.4660 

171.8138 

481.9206 

944.1763 

526.8447 

847.2186 

27.4655 

171.7006 

479.9624 

941.5603 

525.6302 

845 2656 

O'T A r> r ri 

L ( .WOO 

171.657b 

479.1306 

936.7340 

524.9711 

844.2057 

27.4652 

171.6385 

478.7476 

934.2196 

524.5739 

843.5670 

27.4652 

171.6291 

478.5541 

932.8928 

524.3162 

843.1526 

27.4652 

171.6240 

478.4485 

932.1524 

524.1396 

842.8685 

27.4652 

171.6211 

478.3871 

931.7163 

524.0132 

842.6654 

27.4652 

171.6193 

478.3496 

931.4469 

523.9198 

842.5151 

27.4652 

171.6182 

478.3256 

931.2736 

523.8487 

842.4008 

27.4652 

171.6175 

478.3097 

931.1583 

523.7934 

842.3119 

27.4652 

171.6170 

478.2988 

931.0791 

523.7495 

842.2413 

27.4651 

171.6155 

478.2670 

930.8451 

523.4749 

841.7997 


Table II. Convergence of Frequencies for Uniform, p- Version Beam 




fv l? Hz 

fv2i Hz 

fv 3 > Hz 

fvii Hz 

fe i, Hz 

ful- Hz 

6 


27.5953 

270.3419 



577.2137 

928.2170 

10 

i 

27.4734 

173.1390 

906.1908 


525.4401 

844.9599 

14 

2 

27.4652 

172.5710 

490.8033 

2110.4920 

523.5106 

841.8572 

18 

3 

27.4651 

171.6203 

490.0998 

988.0348 

523.4752 

841.8003 

22 

4 

27.4651 

171.6178 

478.4139 

987.4482 

523.4749 

841.7997 

26 

5 

27.4651 

171.6155 

478.4039 

932.4859 

523.4749 

841.7997 

30 

6 

27.4651 

171.6155 

478.2675 

932.4643 

523.4749 

841.7997 

34 

7 

27.4651 

171.6155 

478.2674 

930.8613 

523.4749 

841.7997 

38 

8 

27.4651 

171.6155 

478.2669 

930.8613 

523.4749 

841.7997 

42 

9 

27.4651 

171.6155 

478.2669 

930.8452 

523.4749 

841.7997 

46 

10 

27.4651 

171.6155 

478.2669 

930.8452 

523.4749 

841.7997 

50 

11 

27.4651 

171.6155 

478.2669 

930.8451 

523.4749 

841.7997 

54 

12 

27.4651 

171.6155 

478.2669 

930.8451 

523.4749 

841.7997 

58 

13 

27.4651 

171.6155 

478.2669 

930.8451 

523.4749 

841.7997 

Exact ... 

27.4651 

171.6155 

478.2670 

930.8451 

523.4749 

841.7997 





Table III. Convergence of Frequencies for Modeling Tapered Beam With Uniform Beam Elements 


Ndof 

Nelem 

fvli Hz 

fv2, Hz 

fv3, Hz 

/t>4? Hz 

hi, Hz 

fu l. Hz 

6 

i 

36.7115 




577.2137 

928.2170 

12 

2 

59.6562 

172.5966 

709.7736 

1391.3137 

934.4916 

1255.9186 

18 

3 

71.9805 

187.9347 

455.8576 

1169.0279 

1089.9564 

1348.4876 

24 

4 

77.8991 


441.9430 

893.5114 

1156.4921 

1381.3936 

30 

5 

80.9493 

218.2811 

455.8657 

834.3421 

1189.1758 

1396.3582 

36 

6 

82.6853 



835.1067 


1404.3645 

42 

7 

83.7585 


481.8723 



1409.1399 

48 

8 

84.4657 

234.3112 

489.7455 



1412.2157 

54 

9 

84.9554 


495.3418 



1414.3128 

60 

10 

85.3083 


499.4259 

877.7242 


1415.8068 

66 

11 

85.5707 

239.5654 


883.2469 


1416.9088 

72 

12 

85.7711 

240.5273 


887.4925 

1238.6846 

1417.7450 

78 

13 

85.9276 





1418.3946 

84 

14 

86.0520 

241.8811 



1241.4746 

1418.9092 

90 

15 

86.1526 

242.3676 


895.6784 


1419.3238 

150 

25 

86.5898 

244.4943 

514.7164 


1246.7811 

1421.1123 

300 

50 

86.7753 





1421.8642 

450 

75 

86.8097 

245.5723 

517.4313 



1422.0033 


100 

86.8373 



911.2465 

1249.2228 

1422.1216 


Table IV. Convergence of Frequencies for Tapered, h- Version Beam 


Ndof 

Nelem 

fv l , Hz 

fv 2, Hz 

fv 3i Hz 

fv 4i Hz 

hi, Hz 

/?il i Hz 

6 

i 

86.9704 

267.8185 



1273.8317 

1434.6550 

12 

2 

86.9065 

251.7443 

533.1733 

1203.8627 

1271.8282 

1432.6379 

18 

3 

86.8561 

247.4589 

538.7430 

940.3664 

1263.2856 

1428.3265 

24 

4 

86.8434 

246.2835 

526.0757 

958.9185 

1258.0545 

1425.9792 

30 

5 

86.8396 

245.9340 

521.0733 

934.7022 

1255.1634 

1424.7097 

36 

6 

86.8384 

245.8121 

519.2780 

921.7652 

1253.4584 

1423.9658 

42 

7 

86.8379 

245.7628 

518.5613 

916.4367 

1252.3826 

1423.4975 

48 

8 

86.8376 

245.7400 

518.2391 

914.1021 

1251.6647 

1423.1850 

54 

9 

86.8375 

245.7283 

518.0774 

912.9746 

1251.1633 

1422.9669 

60 

10 

86.8374 

245.7216 

517.9880 

912.3746 

1250.8000 

1422.8087 

66 

11 

86.8374 

245.7176 

517.9345 

912.0268 

1250.5286 

1422.6906 

72 

12 

86.8374 

245.7149 

517.9004 

911.8101 

1250.3208 

1422.6002 

78 

13 

86.8373 

245.7131 

517.8775 

911.6672 

1250.1582 

1422.5294 

84 

14 

86.8373 

245.7119 

517.8615 

911.5687 

1250.0287 

1422.4729 
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Tabic V. Coiivei geiiCe 


/uT Hz 

fv 2, Hz 

86.9704 

267.8185 

86.9555 

249.6871 

86.8675 

248.6826 

86.8417 

246.2815 

86.8376 

245.7577 

86.8374 

245.7151 

86.8374 

245.7151 

86.8373 

245.7127 

86.8373 

245.7101 

86.8372 

245.7087 

86.8372 

245.7080 

86.8372 

245.7078 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 

86.8372 

245.7077 


Frequencies for Tapered, p- Version Beam 


fv 3, Hz 


607.6043 

545.9870 

537.1999 

520.7568 

517.9835 

517.9218 

517.8910 

517.8475 

517.8238 

517.8140 

517.8104 

517.8092 

517.8088 

517.8087 

517.8086 

517.8086 

517.8086 

517.8086 

517.8086 


fv 4) Hz 


1173.2573 

1018.3391 

981.2253 

920.7022 

912.1356 

912.1356 

911.6613 

911.3808 

911.2864 

911.2582 

911.2498 

911.2473 

911.2466 

911.2464 

911.2463 

911.2463 

911.2463 

911.2463 


fe i, Hz 
1273.8317 
1254.6750 
1249.4274 
1249.2122 
1249.2122 
1249.2093 
1249.2074 
1249.2067 
1249.2065 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 

1249.2064 


full Hz 
1434.6550 
1434.2370 
1422.2183 
1422.1270 
1422.1223 
1422.1164 
1422.1149 
1422.1145 
1422.1145 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 

1422.1144 






Physical beam 


p-version model 


h-version model 

Figure 3. Uniform cantilevered beam and representative h- and p-version finite-element models. 
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Number of degrees of freedom 


Figure 4. Percent error versus number of degrees of freedom for first bending frequency of uniform beam. 






Number of degrees of freedom 

Figure 5. Percent error versus number of degrees of freedom for second bending frequency of uniform beam. 



Number of degrees of freedom 

Figure 6. Percent error versus number of degrees of freedom for third bending frequency of uniform beam. 
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Number of degrees of freedom 

Figure 7. Percent error versus number of degrees of freedom for fourth bending frequency of uniform beam. 



Number of degrees of freedom 

Figure 8. Percent error versus number of degrees of freedom for first torsional frequency of uniform beam. 
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p-version 

h-version 


Percent error 



0 6 12 18 24 30 36 42 48 54 60 66 72 78 84 

Number of degrees of freedom 

Figure 9. Percent error versus number of degrees of freedom for first axial frequency of uniform beam. 
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Tapered p-version model 
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Uniform h-version model 


Figure 10. Tapered cantilevered beam and representative h- and p-version finite-element models. 
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Tapered h-version 
Uniform h-version 



Figure 11. Percent error versus number of degrees of freedom for first bending frequency of tapered beam. 



Figure 12. Percent error versus number of degrees of freedom for second bending frequency of tapered beam. 
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Figure 13. Percent error versus number 


Number of degrees of freedom 

of degrees of freedom for third bending frequency of tapered beam. 



Figure 14. Percent 


error versus number of degrees of freedom for fourth bending frequency of tapered beam. 
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Figure 15. Percent error versus number of degrees of freedom for first torsional frequency of tapered beam. 



Figure 16. Percent error versus number of degrees of freedom for first axial frequency of tapered beam. 
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